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Abstract 

We present a new derivation of the Onsager-Joyce-Montgomery 
(OJM) equilibrium statistical theory for point vortices on the plane, 
using the Bogoliubov-Feynman inequality for the free energy, Gibbs 
entropy function and Landau's approximation. This formulation links 
the heuristic OJM theory to the modern variational mean field the- 
ories. Landau's approximation is the physical counterpart of a large 
deviation result, which states that the maximum entropy state does 
not only have maximal probability measure but overwhelmingly large 
measure relative to other macrostates. 
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1 Introduction 



The so-called mean field equation of the Onsager- Joyce-Montgomery (OJM) 
theory^ '^'^ for the equilibrium vorticity distributions of a two-dimensional 
inviscid and incompressible fluid (and guiding center line charge model in 
plasmas) were derived more than twenty years ago but the volume of re- 
search concerning these equations continues unabated. Recently, the mean 
field thermodynamic limit for the planar vortex model has been shown to 
exist by Caglioti et al and independently by Kiessling^, and the existence 
of negative temperature states in this model was demonstrated by Eyink 
and Spohn ^. Furthermore, these authors established the maximum entropy 
principle^ for this theory. Onsager^ predicted that the negative temperature 
equilibrium states exhibit a clustering of the point vortices into coherent 
structures. Exact solutions of the OJM mean field equations obtained by 
Chen et al ^ confirm the existence of such coherent states in inviscid turbu- 
lent flows. 

A OJM type equilibrium statistical theory for nearly parallel thin vortex 
filaments was very recently put forth by Lions and Majda^. They established 
the mean field thermodynamic limit and the maximum entropy principle for 
their problem. The mean field equation that Lions and Majda obtained, 
is radically different from other mean field theories for vortex dynamics in 
the sense that their equations are time-dependent. Nonetheless, in a special 
singular limit of perfectly parallel vortex filaments, their mean field equation 
reduces to the sinh-Poisson mean field equation of the planar point vortex 
problem. 

Concerning the relationship between the OJM and other statistical theo- 
ries, Majda and Holen^° showed that both the OJM theory and Kraichnan's 
truncated spectral theory^^ are statistically sharp with respect to the recent 
infinite constraints theory of Miller^^'^^ and Robert^^, i.e., the few constraints 
and infinite constraints theory agree for low energies. Chorin^^ has earlier 
indicated that a few constraints is sufficient for a reasonable theory of sta- 
tistical equilibrium in these models. He also demonstrated numerically that 
the Miller-Robert theory is valid only for moderate temperatures and has no 
Kosterlitz-Thouless phase-transition except at zero temperature^^. Turking- 
ton and Whittaker^ gave a numerical algorithm for the Miller- Robert theory. 
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Recently, Turkington^^ critiqued the Miller-Robert theory and presented a 
few constraints equilibrium statistical theory for 2-D coherent structures, 
where the equality constraints of the Miller-Robert theory are replaced by 
inequalities. Turkington^*^ argued that the finite ultra-violet cutoff implicit 
in the Miller-Robert theory does not reflect the true inviscid vortex dynamics 
as represented in the two dimensional Euler equations. In their preprint 
Boucher, Ellis, and Turkington presented rigorous large deviation results and 
proved the existence of the mean field thermodynamic limit for Turkington's 
theory. 

Yet, several issues remain unresolved — in spite of the above recent works, 
it is still not completely clear where the OJM theory stands in relation to 
other equilibrium statistical theories of two dimensional turbulence^^'^^'^^'^'^ 
especially with regard to their relative efficacy in modelling real physical 
problems in fluid mechanics. Moreover, it is still not clear how the OJM 
theory is related to the physicists' standard variational formulations of mean 
field theory such as in Chapter 6 and 7 of the text^''. It is the main aim of this 
paper to derive equation (0) in a new way, thereby demonstrating that the 
OJM theory is indeed a mean field theory in the sense of Weiss, Peierls and 
Feynman^° . Specifically, we will prove that equation (|1|) is the limit, as the 
number of particles iV ^ oo, of mean field equations that are derived using 
the Bogoliubov-Feyman inequality^^ (Q) , the Gibbs entropy function and 
Landau's approximation^^. Previous derivations of the OJM theory^'^'^ are 
based on Boltzmann's entropy function instead of Gibbs' entropy function. 
We note here the significant point that our result holds for all initial vorticity 
distributions qo{x). 

A OJM-type theory for point vortex systems on a rotating sphere^^ has 
recently been obtained by Lim^^, and the corresponding mean field thermo- 
dynamic limit was established by methods similar to those of Kiessling^. 

2 Mean field theory 

All the above statistical theories for two dimensional turbulence in inviscid 
fluids have a common thread in the maximum entropy principle of informa- 
tion theory''. They lead to nonlinear elliptic equations, the so-called mean 
field equations which have the generic form 

= F(^,7) 
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for the stream function \1/ of the mean field locally averaged vorticity dis- 
tributions q{x) = — A^, and some parameters 7. Indeed it can be shown 
that q{x) are exact stationary solutions of the Euler equations. The mean 
field equation of the OJM theory for a single component vortex gas is the 
nonlinear elliptic equation 

A* = A;e-^*. (1) 

Heuristically, this equation has the correct form for a mean field theory in 
the sense that the precise interactions between particles has been replaced 
by an approximate energy expression where the particles interact with some 
mean field. 

The main tool for the work in this section is the Bogoliubov-Feynman 
inequality for the free energy in the form proven by Feynman^^: 

F<Fo+{H,)^ = F,ar, (2) 

where Fq is the free energy based on the approximate Hamiltonian Ho, where 
H = Hq + Hi, and the second term on the right is the thermal average of 
the remainder Hi in the canonical measure based on Hq. In other words, 

_ JcdV (•) eM-PHo} ^ 

where the partition function Zq = J^, dV exp{— /3/7o} is based on the approx- 
imate Hamiltonian Hq. In the standard mean field theory^", one choose Hq so 
that Zq can be evaluated easily, and yet it must have certain basic physical 
properties to make the theory physically relevant. For negative temperatures, 
a derivation similar to Feynman's yields 

F>Fq + {Hi), ^ F-,. (4) 

For the OJM theory of point vortices in the plane where the vortices are 
identical and have vorticity or charge A, we shall choose the coarse-grained 
approximate Hamiltonian based on the division of the physical domain Q C 
i?^ of area A into M equal boxes. 



M M 

-T^Y^Yl niTijX^ log - f° 

^ i=l j^i 
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Here ni denotes the number of vortices in box Bi which has area K^, and 
is the location of the center of Bi. Furthermore, the total number of particles 
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is N, i.e., 



M 



Thus, 5° are no longer dependent on time but depend instead on the lat- 
tice which is implicit in the coarse-grained Hamiltonian Hq. Since the full 
Hamiltonian of the point vortex model is 
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N N 



the remainder Hi is given by 

Hi — H — Ho 

1 M I Hi rii 



M M 

EE 

1=1 i'^l 
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1=1 \j=i fc^j 

E E log \xj - ffc I I - n^ni/ log |f° - | 

-I M / Ui Tli 

-oE EEA^iogl4-4l 
^ i=i \j=i k^j 



1 



Af M 



-^EE 

^ j=i iVi 



^ ^ log I (f° - f«0 + (4 - 4) I - n^n,\' log - 5?°, | 

\j=l k=l J 

where we have used 

-> -» / -0 -0\ if-'/ 



in the first sum, and 



Xj — Xk — — x^i) -\- yXj — x^) 



in the second sum. In other words the primed vector x'j denotes the difference 
Xj — x^, which represents the vector from the center of Bi to vortex j in it. 
From 



one derives 

iogi(x°-4)+(4-4)i - iogi(f°-4)i+iog(i+ [||:||| ) 

using the complex form of the logarithm, a step which we omit here for the 
sake of brevity. Substituting this last inequality back into Hi and keeping 
only terms of order \{x'j — x'^)\/\{^ — x^,)\, we get 

M [ Ui Ui \ M M ( rii n^i i/-*/ _ -j/ x 

^ i=\\j=\k^j ) ^ j=i j/^i W-i^i -i^i', 

where the first group of terms represent the intra-box interaction energy, and 
the second group represents 0{K) terms in the inter-box interaction energy. 

To begin, we compute the partition function Zq which is based on the 
approximate Hamiltonian i^o, as follows 



Zo = ^iy(s)/i'^exp(-/?i7o) 

s 

= Y: W{s)h''' e^pl^EJ: log - 41 , 

s V 31^i J 

where the summation is taken over all coarse-grained state (or macro-state) 
s — (ni, ...^um) such that Z^^i/ij = and the degeneracy of the macro- 
state s is given by 

N\ 

W{s) ^ — -. 

ni\...nM\ 

The probability distribution Po{s) for a macrostate s, is now defined by: 

P„M = j J^^^i^A (5) 

W{s)h^^exp{-pHo) 



where D^{s) denotes the part of phase space which is occupied by the 
microstates associated with s. Finally, we compute the free energy Fq as 
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follows: 



Fq = --logZo 



1 

1 



M M 



N\ogh^ + \og\J2W{s)exp ^^^nin,A2log|x°-:^ 



In the next step in this procedure, we will compute {Hi)q up to 0{h) : 



(^i)o 



HAx) — — ^ — ^^^^^ dxi...dxN 



M ( n-i Hi 



hi M Tii n^i 

4E E E E A 
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1 M I rii Ui 



X i- 



exp {-f3Ho{x)) 
exp (-/3i7o(^)) 



dxi...dx 



N 



dxi...dxN +0{h). 



where x = (xi, xa?) is a microstate. Substituting the above expressions for 
Zo, Fq into (H), the free energy upper bound when the temperature P > 0, 
is given by 

Fvar = -fo + {Hi)q 

-iiVlog - i log 1^ W{s) exp (^f E E n,n,A2 log - | 



+ 



|E ( E E A^loglf, 



«=1 \j='i-k^j 



exp(-/3Ho(i^)) 
^0 



dxi...dxN 



similarly for negative temperatures /? < 0, the free energy lower bound is 
given by 

Kar = ^0 + (^l)o 



_ 1 AT log - ^ log E W^(^) exp f E E 

I s \ i=^j¥'i 



niUjy log |x° — x*- ' 



+ 



M / Hi rii 

-|E E E A^loglx;.-^ 



Whether the temperture is positive or negative, the integral in the above 
equation can be evaluated as follows: 



M ( TLi Ui 

E E A'log|f;-4 

(M M 
f x;E/^i%^^i°gi^-^i 
i=l jyti 



X 



Zo 



(6) 



2Zo 



s 



/MM 

exp f E E niTijX log |a?^ — 5° 



X 



M r 

^ / 

i=l ^ 



2 ^ '"vy 

E E A'log|f' -41 

i=i fc^i 



2^0 



E 



Mr ( rii -Hi \ 

X E / c?'''^' E E A'log|f' -fj^l 



> . 



An application of the Mean Value Theorem implies that the integrand in 
each of the integrals in the last hne above can be replaced as follows: 



EEA^iogi4-4l 



y^Uiijii - 1) 



L(ni), 



where L{ni) is a large negative constant of the order of 

h 



log 14 - 41 - log 



Therefore, the term 



M . ( Ui Ui 



(7) 



is given by 



1=1 



Substituting this expression in (13), we get 



1 



1 M I rii rii 

4e EEA^iogl4 



exp {-PHo{x)) 



2Zf 




dxi...dxN 



L{ni)K 



2ni 



'0 S \^ J 4 = 1 

Thus, both F~^^. and F^ar up to order 0{h) is given by: 



^var ~ Pvar 



1 
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/5 



M M 



N log + \og{Y^ W{s) exp I I 5] 5] n,n, log jx" - 



2 \ M 



Y.X^n,{ni-l)L{n,)h^''^)+0{h) 



i=l 



2.1 Bounds 

The strategy of variational mean field theories^° at this point is to intro- 
duce some parameters 7 into the expressions for F^ar (respectively F'^^^) and 
minimize F^ar{l} (resp. maximize F~^^{'~^}) with respect to 7. This yields 
the best possible approximation to the full free energy F within the class 
determined by the chosen reduced Hamiltonian Hq. Since 

we now have the best approximation for Z. In our derivation of the mean 
field theory for the point vortex model, these parameters are given by the 
macrostate s = {rii, ...um} which are occupation numbers of the boxes Bi in 
the statistical coarse-graining procedure. We will now use the facts that the 
energy H of the system, the entropy Sq and the partition function Zq are 
all highly focussed at the equilibrium (most probable) macrostate s* in the 
mean field limit of large N. This is essentially the Landau approximation in 
the Ginsburg-Landau theory for phase transitions^^. In particular, we will 
use the following consequence of the Landau approximation 



W{s*) 


/ M 

n 


A- 






( M 








n 






s 
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to analyse the expressions (j^). Thus F^ar and F^^^. become 



F~ = F ~ 

var yar — 



(3 



M M 



N log + log { Wis*) exp ^ ^ ^ mnjX^ log 



1^ _ -^J0| 



=1 

~4 ' 



M 



^ A^ni(ni - l)L{ni) 



for all /5 > and for a range of values > /5 > where is given by the 
theory for the existence of the mean field thermodynamic limit 

In order to minimize F^ar and maximize F^^^ with respect to s = 
{rii, ...,nm} under the constraint 



M 



i=l 



we augment F^aris) and by adding the auxiliary term with Lagrange 

multiplier a to obtain 



M 



Fr,,^^ [ Sj 



+ a{Y,ni-N), 

1=1 

M 



i=l 



Then taking the gradient V sF^ar{.s) = (resp. VgF^ari^) = 0) yields 



M 



/3 



{logUi + 1) + ^ A\log 



^ _ ^1 

•^i -^j I 



ni(jk-l) , 2n,-l ni(jii-l) 
1 L [ni) + — - — Lim) + Lim 



+ a = (10) 



for i = 1,...,M. We note that VsF^aris) (resp. VsF^^^{s)) differs from 
VsFo(s) by the derivatives of the term 



E 

j=i 



L{ni), 



which is associated with the self energies of each of the M boxes Bi in the 
coarse-graining procedure. In working with F^aris*) (resp. F~^^{s)) instead 
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of Fq, we have put back these self energy terms. Minimizing F^ar (resp. max- 
imizing F~^^ in the case of negative temperatures) yields a lower bound for 
Fvar (resp. an upper bound for F~^^) which provides the best approximation 
for the free energy F of the point vortex system^° within the ansatz of as- 
suming Ho to be the approximate (mean field) Hamiltonian. Solving for Ui 
in (p!OD yields the occupation numbers 



M 



Hi 



e"" exp I -/5 I E log |x° - + 



riiiui - 1) , 



L'irii) + 



2nf - 1 



-Lirii 



Next we turn to the analysis of the behaviour of this expression as A^, M — >■ 
oo while the vortex strength is scaled by j^, i.e., A = — it is easy to see 
that the following limits are valid: 



1 



N 

—a 



riiirii - 1) , 
L (rii 



2n? 



-L(nj) 



(12) 



M 



^log|a;, -x^\ 



Thus the mean field equations in the planar point vortex theory is given by: 



C(x) = = rfexp (-/3(E°(x) + E\x)) 



(13) 



for all temperatures, which differs from that in the OJM theory by the self 
energy term E^{x). 

We will now analyse the self energy expression in some detail. As we 
let N, M —>■ oo the vortex strengths have to be scaled by in order that 
the mean field (non-extensive) thermodynamic limit exists, as was shown by 
Caglioti et al^, Kiessling^ and Eyink and Spohn^. Since < Mmaxj(?7,j) 
and h"^ = A/M, we have the bound 



w ^ I, 1, 



max,- Tij 
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= J log A - ^ log M - J log[max(ni)] 
z / / « 

= -logA - -log[Mmax(ni)] 
> -logv4 logA^, 



which imphes that for each i 

\L{n, 

The self energy term in (|T3|) scales as follows 



< ^logAT-^logA 



(14) 



\E\x)\ 


< 




L' 




1 + ^1^(^01 




< 




L' 




Hi maxj rii 

1+ iV 




< 


9 

nf 


L' 







where we have used the fact that 
as N ^ oo because 



7V2 



\L'in.)\ 



< 



N 

1 

2N 



< 1. The first term tends to zero 



1 

2n,; 



2m 



0. 



By (|Tip, the second term tends to zero as 
\L(n, 



rii 



N 



< 



N 



1 1 
--logA + -logAr 



because the number ~ of particles in box Bi stays about the same as 
the total number of particles and the number M of equal boxes in the 
statistical coarse-graining procedure both tend to oo. We have shown that 

Theorem: The mean field equation (13) of the Bogoliubov-Feynman- 
Landau mean field theory for point vortex dynamics tends to the mean field 
equation (|l]) of the OJM theory in the mean field thermodynamic (non- 
extensive) limit of infinite particles N oo, and infinite number of boxes 
M — i> oo in the coarse-graining procedure in the definition of the approximate 
Hamiltonian Hq. 
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3 Concluding remarks 



The following discussion gives a brief derivation of another formulation for 
F^ar and F'^^^. By definition, the free energy based on the Hamiltonian Hq 
is given by 

-^0 = (-f^o)o ^ TSo, 

where the expectation operator is defined by (|^), and 5*0 is the Gibbs 
entropy function based on i^o, i-e., 

s 

with 

p , , _ iy(5)/i^^exp(-/3i7o(3)) 
^o\S) = . 

^0 

Thus, the upper bound F^ar (resp. lower bound F~^^) for the free energy F 

^var = Pvar = {Hq + Hi)^^ — TSq. (15) 

is exactly equal to the free energy based on the full Hamiltonian H, but using 
the probability distribution Pq. 

To summarize, we have derived the OJM theory by a new method which 
is based on the Bogoliubov-Feynman inequality, the Gibbs entropy function 
and Landau's approximation, and showed that it is a mean field theory. It is 
remarkable that no additional conditions on the initial vorticity distributions 
was needed to prove this result. In the case of a two component vortex gas, 
our procedure gives the well-known sinh-Poisson equation^. The analogue 
of the above theorem for point vortices on a rotating sphere is presented in 
Lim^^. We note that previous derivations of the OJM theory^'^'^ are based 
on Boltzmann's entropy function instead of Gibbs' entropy function, which 
is the fore-runner of the information-theoretic entropy function. 

In another paper^^, this author will return to the issue of the indetermi- 
nacy of the OJM mean field equations which was raised by Turkington 



This issue concerns the fact that a given continuous vorticity distribution can 
be represented in a number of different ways by clouds of point vortices — 
for example, one could use two species of vortices with equal but opposite 
circulations, or one could just as well choose an approximation based on 
three different species of vortices. The mean field equations ensuing from 
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these distinct representations of the original continuous vorticity distribu- 
tions must necessarily differ; this can be demonstrated equally well within 
the traditional formulation and the current derivation of the OJM theory. 
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